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Abstract 


"  In  an  earlier  paper  (Shaub-  numerical  solutions  of  the  Landau 
equation  were  obtained  which  show  that  the  tail  of  the  electron  velocity 
distribution  function  differs  substantially  from  a  local  Maxwellian  distribu¬ 
tion  in  the  solar  transition  region  and  upper  chromosphere.  In  this  paper,  I 

show  that  a  linearized  version  of  the  BGK  model  kinetic  equation,  with 

1 

collision  frequency  proportional  to  can  solved  analytically  for  the 
tail  of  the  distribution  function  in  an  atmosphere  with  prescribed  temperature 
and  density  profiles.  Results  for  the  angle-averaged  distribution  so  obtained 
are  shown  to  be  in  reasonably  good  agreement  with  earlier  numerical  results. 
Accurate,  easily  evaluated  approximations  for  the  tail  of  the  distribution 
function  are  derived  from  the  exact  formulas.  These  show  that  both  the  zeroth 
and  first  angular  moments  of  the  distribution  function  are  nearly  power  laws 
over  a  wide  velocity  range  in  the  low  transition  region.  I  also  show  that  the 
heat  flux  into  lower  temperature  region  is  carried  by  suprathermal  electrons 
with  velocities  well  above  local  thermal  velocities.  The  formulas  given  here 
should  be  useful  in  the  calculation  of  electron-ion  inelastic  collision  rates 
under  conditions  in  which  the  local  Maxwellian  approximation  is  invalid. 
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I.  Introduction 


In  Shoub  (1983;  hereafter  referred  to  as  Paper  I)  numerical  solutions  of 
the  Landau  equation  were  obtained  for  the  electron  velocity  distribution 
function  (EVDF)  in  idealized  models  of  the  solar  transition  region.  These 
results  showed  that  the  EVDF  forms  a  pronounced,  anisotropic,  high-velocity 
tail  in  the  low  transition  region  (!  <  3  *  10^  R)  and  upper  chromosphere  as 
the  result  of  the  streaming  of  fast  electrons  down  the  temperature  gradient. 
It  was  also  shown  that,  in  consequence,  inelastic  electron-ion  collision  rates 
are  signficantly  enhanced  in  these  regions  relative  to  their  values  obtained 
using  a  local-Maxwellian  energy  spectrum. 

The  results  reported  in  Paper  I  were  obtained  through  involved  numerical 
calculations  whose  nature  and  expense  hinder  the  development  of  physical 
intuition  concerning  the  underlying  transport  processes.  In  this 
complementary  paper  I  therefore  address  the  same  problem  as  considered  in 
Paper  I,  but  here  use  a  Bhatnagar-Gross-Krook  ^(1954)  model  for  the  collision 
term.  Moreover,  I  take  advantage  of  my  earlier  finding  that  only  the  tail  of 
the  EVDF  differs  significantly  from  a  local  Maxwellian  under  transition  region 
conditions  in  order  to  linearize  the  BGK  collision  term  and  thereby  obtain  an 
analytically  soluble  kinetic  equation.  In  this  approximation  one  can  obtain 
an  expression  for  the  EVDF  in  an  atmosphere  with  a  prescribed  temperature  and 
density  profile.  Results  for  the  angle-averaged  distribution  function  so 
obtained  are  shown  to  be  in  reasonably  good  agreement  with  those  found  earlier 
using  the  Landau  equation. 

Although  the  analytical  solutions  I  obtain  are  relatively  complicated 
[see  equations  (7a)  and  (7b)],  they  can  be  accurately  approximated  at  supra- 
thermal  velocities  by  much  simpler  expressions  which  provide  insight  into  the 
underlying  physical  processes.  These  formulas  will  hopefully  be  useful  to 
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workers  wishing  to  evaluate  electron-ion  inelastic  collision  rates  under 
conditions  in  which  electrons  do  not  have  a  local  Maxwellian  energy  spectrum. 

Two  new  results  which  emerge  from  the  present  work  are  a)  that  the 
angle-averaged  distribution  function  is  nearly  a  power  law  over  a  wide 
velocity  range  in  the  lower  temperature  regions  of  the  atmosphere,  and  b) 
that  the  conductive  energy  flux  into  the  lower  temperature  regions  is  carried 
predominantly  by  suprathermal  electrons  with  velocities  3  to  25  times  local 
thermal  velocities.  The  latter  result  shows  clearly  that  classical  transport 
theory,  which  predicts  that  the  heat  flux  is  carried  by  electrons  with 
velocities  two  to  three  times  local  thermal  velocities,  is  invalid  under 
conditions  met  in  the  solar  transition  region. 
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II.  Formulation 


Consider  an  inhomogeneous  but  isobaric  slab  of  fully  ionized  hydrogen  in 
which  the  protons  are  idealized  as  being  infinitely  massive  and  at  rest.  The 
slab  has  thickness  L,  and  any  magnetic  field  present  is  assumed  to  be  parallel 
to  the  gradient  direction,  which  I  take  to  be  the  z-axis.  Let  f(u  ,v,z) 
denote  the  EVDF,  with  v  the  electron  velocity,  and  y  *  v_*e z/ 1 v.1 .  The  form  of 
the  BGK  equation  I  shall  use  is  then 

df  * 

UV  —  *  y(v,z)  {f  (y  ,v,z)  -  f (y  ,v,z)>  ,  *  (In) 


where 


f*(M,v,z)  *  n'(z)(m/27rkT'(z))3^2  exp[-m(v-w'ez)2/2kT'(z)]  (lb) 


is  a  local-Maxwellian  distribution  with  parameters  n',w',  and  T'  and  (cf., 
Rawls,  Chu,  and  Hinton  1975),  and 


v (v, z) 


16tt e^nlnA  ,  .  2  , 

(erf(£)  -  y=-  £exp(-£z) } 


m2v3 


(lc) 


with  E,  ■  v/vth  and  v*h  *  2kT/m.  Here  m  is  the  electron  mass,  InA  is  the 
coulomb  logarithm,  and  F  is  a  constant  to  be  chosen  later.  When  r  ■  1,  v 
is  the  slowing  down  time  of  a  test  electron  in  a  fully  ionized  hydrogen  plasma 
(Krall  and  Trivelpiece  1973). 

Several  remarks  concerning  equation  (la)  are  in  order.  First,  note  that 
although  at  first  sight  (la)  appears  to  be  a  linear  equation,  it  is  not.  In 
order  to  conserve  particles,  momentum,  and  energy,  the  parameters  n',  T  and 
w'  appearing  in  f*  must  be  chosen  to  satisfy  the  constraints 
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particles 
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Che  corresponding  thermodynamic  variables  which,  in  turn,  may  be  estimated 
using  macroscopic  (i.e.,  classical)  arguments.  I  therefore  choose 


T'(z)7/2 

-  I1'1  +■  (t£/2  -  T7/2)  z/L  , 

(3a) 

n'T' 

*  constant  , 

(3b) 

* 

w 

-  o  , 

(3c) 

corresponding  to  an  isobaric  slab  of  thickness  L  in  which  energy  is 
transferred  solely  by  (classical)  thermal  conduction.  (To  arrive  at  (3a)  I 
have  treated  In  A  as  a  constant.]  Here  and  Tc  are  temperatures 
characterizing  the  incoming  EVDF  at  the  upper  and  lower  boundaries, 
respectively.  With  n',  T'  and  w'  regarded  as  known  functions  of  position, 
equation  (la)  reduces  to  a  linear,  ordinary  differential  equation  in  which 
U  and  v  appear  as  parameters.  This  represents  a  significant  mathematical 
simplification  relative  to  the  Landau  equation,  where  U  and  v  enter  as 
independent  variables.  Note,  however,  that  because  the  above  constraints  are 

I 

violated  in  this  approximation,  one  cannot  obtain  a  kinetic-theoretic  predic¬ 
tion  of  the  electron  temperature  profile.  One  must  solve  the  nonlinear 
problem  for  this. 

In  order  to  retain  the  advantage  mentioned  above,  I  have  omitted  a 
thermal  electric  field  term  from  the  lefthand  side  of  (la).  This  is  not  a 
serious  omission  for  my  present  purposes  because,  as  argued  in  detail  in 
Paper  I,  this  term  is  unimportant  relative  to  the  streaming  term  at 
suprathermal  velocities. 

It  is  useful  to  note  that  by  analogy  with  the  equation  describing  photon 
transport  in  a  medium  in  which  the  matter  is  in  a  state  of  local  thermodynamic 
equilibrium,  the  linearized  BGK  equation  can  be  interpreted  in  the  following 
simple  way.  Electrons  travel  linear  trajectories  at  constant  velocity,  with 
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Che  pathlength  of  a  given  flighC  being  a  random  variable  whose  probability 
distribution  is  fixed  by  Che  electron's  speed  and  collision  frequency  with 
background  electrons  and  protons.  A  given  flight  is  terminated  in  a  single 


thermalizing  event  in  which  the  electron  is  "absorbed"  and  then  "re-emitted" 
with  a  new  random  velocity  whose  probability  distribution  is  f* .  This,  of 
course,  is  not  how  electrons  behave.  Nevertheless,  proper  choice  of  the 
collision  frequency  v  allows  many  aspects  of  more  complicated  descriptions 
(e.g.,  via  the  Landau  equation)  to  be  modeled  accurately. 

Since  my  interest  lies  in  the  high-velocity  tail  of  the  distribution,  one 
final  simplification  I  shall  invoke  is  to  use  the  large-velocity  form  of  the 
collision  frequency  given  in  equation  (lb).  It  is  not  necessary  to  do  this, 
but  the  resulting  formulas  are  more  transparent  if  I  do.  Further,  I  shall 
omit  the  primes  from  n',  T'  and  w',  because  I  no  longer  distinguish  between 
these  variables  and  their  thermodynamic  counterparts. 
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Ill .  Analysis 


a)  Formal  Solution 

Let  us  begin  by  introducing  x  =  InT  as  the  independent  spatial 
variable  and  use  (1c)  and  (2b)  to  rewrite  (la)  in  the  form' 


tide 


df 

dx 


f*  -  f 


(4a) 


where 


a(T)  -X(T)  ;  \  (T)  -  (kT)2/4r^e4nln:\ 

dz 

5(T)  »  v/vth(T)  ;  v2h  *  2kT/m  . 

It  follows  from  (3a),  (3b)  and  (4b)  that 


(4b) 

(4c) 


a(T) 


2  M1/2 

7  L  \T  / 


(4d) 


Equation  (4a)  is  to  be  solved  subject  to  the  boundary  conditions 

f(z*0)  ■  f*(v;Tc,nc)  *  nc  (m/2TrkXc)2^2  exp[-(mv2/2kTc)] ;  >.  0  (5a) 
f(z-L)  -  f*(v;Th,nh)  -  nh  (m/2rrkTh)3/2  exp(-(mv2/2kTh) ] ;  <0  (5b) 


with  "  ncTc.  These  are  the  same  boundary  conditions  as  used  in  Paper  I. 

The  choices  T ^  «  2  x  10®  K,  Tc  ■  8.1  x  103  K,  n^  »  3  x  10®  cm  3  and 

Q 

L  *  5  x  10  cm  will  hereafter  be  referred  to  as  the  standard  model. 

The  solution  of  equation  (4a)  subject  to  conditions  (5a,b)  may  be  written 


as 


8 


1 


f<U,v;T)«  £*(Tc)exp[-pv(T,Tc)/ U] 


T 


•/ 


+  /  f  (T') 


exp  [  -p  (T,T')  AJ  dT' 

- 7 - ;  o<jj<i 

ya  (T')S(T')4  T' 


(6a) 


and 


f(U',v,T) 


f*(Th)exp[-pv(T,Th)/u'] 


T 


exp[-pv(T,r)Xi  ']  dr 

- - - ;0^j'<1  (6b) 

ua(I')((I')4  r 


where  u 


-y  and 


PV(T,T') 


_ 2 

SqlCt'X  (T')4 


111 


(6c) 


The  quantity  exp[-pv(T,r)/ 1  u I )  is  the  probability  that  a  (u  ,v)  electron 
"emitted"  at  T'  will  survive  a  flight  to  T.  Note  that  for  fixed  v,  py  is 
symmetric  in  T  and  T' ,  since  ct?4  is  proportional  to  T  An  electron  of 

given  velocity  is  thus  equally  likely  to  go,from  T  to  T'  as  from  T'  to  T. 
However,  it  is  much  more  likely  that  an  electron  "emitted"  with  a  given  value 
of  5  will  survive  a  flight  to  lower  temperatures  than  one  to  higher  tempera¬ 
tures.  This  is  an  important  point,  because  at  any  given  location  electrons 
are  available  only  over  a  limited  range  of  £> 

It  will  be  helpful  later  on  to  notice  that  equation  (4a)  implies  that 
the  quantity  ct?4  may  be  interpreted  as  the  electron  mean  free  path,  measured 
in  units  of  InT. 

Equations  (6a, b)  may  readily  be  integrated  over  y  to  obtain  the  zero'th 
and  first  angular  moments.  When  written  in  dimensionless  form,  the  results 
are : 
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1/v i{(T/Th)5/2  exp(-C2  T/Th)  E2[pv(T,Th)] 
+  (T/Tc)5/2  exp(-' 2  T/Tc)  E2[ pv(T>Tc) ] 


T 

,V‘/ 


exp(-- 2  T/T')  E1[pv(T,T')]  dT'/T'}  , 


=  7,7  ((T/Tc)5/2  exp(-52  T/Tc)  E3  [pv(T,Tc>  ] 


-  (T/Th)5/2  exp(-^2  T/Th)  E3[pv(T,Th) 


♦  (ctC4)-1  I  exp(-£z  T/T')  E2tpv(T,T')]  dT'/T' 


-h 

-  (a C4)"1/  exp(-£2  T/T')  E2(pv<T,T')]  dT'/T'}  . 


Here  <J>  -  (2nv  J/n)  f  ,  i>°  -  1/2 


1  1 
jt>  dV  ,4>X  -  3/2  Jv4>  du. 


ao 

Eq(x)  *  t-n  exp(-xt)dt,  n  *  1,2,3,...,  are  exponential  integral 
functions . 

The  right  side  of  equation  (7a),  evaluated  with  r  *  1  and  parameter 
values  equal  to  those  of  the  standard  model,  is  compared  to  corresponding 
Fokker-Planck  (FP)  results  in  Figure  1.  It  is  seen  that  although  BGK  theory 
(with  r  *  1)  predicts  somewhat  larger  deviations  from  a  local  Maxwellian  than 
does  FP  theory,  the  results  are  qualitatively  similar,  at  least  over  the 
velocity  range  in  which  the  FP  results  have  been  calculated.  I  find  that 
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choosing  r  ■  2.2  gives  quite  good  agreement  in  the  tail  of  the  distribution. 
Note  that  since  Figure  1  is  a  log-log  plot,  the  nearly  linear  behavior  of 
curves  at  intermediate  velocities  implies  that  3®  is  approximately  a  power  law 
over  this  velocity  range.  This  result  is  derived  below. 

Figure  2  shows  the  dimensional,  angle-averaged  distribution 
f°  ■  (27Tvth^  n't®  at  several  locations  with  corresponding  local-Maxwellians 
shown  as  dashed  lines  for  comparison.  The  relatively  weak  spatial  (tempera¬ 
ture)  dependence  of  the  high-velocity  tails  should  be  noted. 

A  Chapman-Enskog  analysis  of  equation  (4a)  yields  the  result  f(y,v,T)  ■ 
f*  +  f^j  +  0(,2),  where 

f] [:l  df*/dx  =  ctC  4(|  ~  .  (8a) 

Thus,  the  classical  BGK  heat  flux  is 

qcl  ”  1/2  m  (2  )  2/3  j  v^  f^  dv  (8b) 

0 

00 

-  4/3  tt“1/2  (nkT)vth  a(T)  J  59(J  -  52)  exp(-£2)  d£  .  (8c) 

o 

The  classical  and  exact  heat-flux  integrand  [as  obtained  from  equation  (7b)] 
are  compared  at  several  locations  in  the  standard  model  in  Figure  3.  From 
these  results  it  is  clear  that  the  conductive  energy  flux  into  the  lower 
temperature  regions  is  carried  predominantly  by  suprathermal  electrons. 

The  heat  flux  obtained  from  (7b)  and  normalized  to  its  classical  value  is 
shown  as  a  function  of  temperature  in  Figure  4.  Results  for  several  pressures 
are  shown;  all  other  parameters  have  their  standard  values.  Since  the 
temperature  profile  I  have  assumed  is  derived  from  the  condition  that  the 
classical  flux  be  constant,  the  spatial  variation  of  the  ratio  shown  in 
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Figure  4  is  due  entirely  to  the  variation  of  the  “Actual"  heat  flux.  That 
this  "actual"  heat  flux  exhibits  a  spatial  variation  is  of  course  incorrect, 
and  is  due  to  violation  of  energy  conservation  (equation  (2c))  in  the  present 
linear  theory.  That  the  calculated  flux  is  not  constant  implies  that  the 
correct  temperature  profile  should  be  different  from  the  T"  law  assumed.  In 
fact,  the  spatial  variation  of  the  flux  shown  in  Figure  4  implies  that  the  gas 
is  cooled  at  high  temperatures  and  is  heated  at  low  temperatures.  Thus,  the 
correct  temperature  profile  will  be  steeper  than  the  classical  profile  at  high 
temperatures  and  less  steep  at  low  temperatures.  Unfortunately,  determination 
of  the  correct  "constant-flux"  profile  requires  solution  of  the  nonlinear 
problem  and  for  this  reason  is  deferred  to  a  later  publication. 


b.  Approximate  Solutions 

The  expressions  for  the  distribution  function  given  above  (i.e., 
equations  (6a, b)  and  (7a, b))  are  difficult  to  evaluate  and  to  interpret.  In 
this  section  I  derive  useful  approximations  to  these  expressions. 

Consider  the  integral  in  equation  (6b).  Using  (2b),  (3a, c)  and  (4c, d), 
it  can  be  rewritten  in  the  form 


J(y',v,T) 


CM(T)(y'a£4) 


TH 

-/ 


i  dT' 

exp[-£2  h(T')]  — 


(9a) 


where  u'  -  -y ,  CM(T)  -  n(T)(m/2  TTkT)'3/2  (Note:  CM(T)  a  T_5/2)  ,  and 


h(T') 


(9b) 


Here,  as  elsewhere,  temperature-dependent  quantities  written  without  arguments 
are  to  be  evaluated  at  T.  The  quantity  hCD  has  a  minimum  at  T  -  T*,  where 
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7/2 

n  • 


(10) 


Note  that  since  both  &(T)  and  £(T)  are  proportional  to  T-^2,  it  follows  from 
(10)  that  T*  is  independent  of  T;  i.e.,  T*  ■  T*<LT,v).  This  point  will  prove 
significant  later  on.  For  large  £  the  integrand  in  (9a)  is  sharply  peaked 
about  T'  *  T  ;  the  integral  may  therefore  be  accurately  approximated  by 
expanding  h(T')  in  a  Taylor  series  about  T'  *  T*,  i.e.,  h(T')  Z  h(T*)  + 
l/2h"(T*)  (T'  -  T*)2,  provided,  of  course,  that  T  <_  T*  <_  T^.  Substituting 
this  expansion  into  (9a)  and  making  use  of  the  relations 


h(T*) 


/T  \ 

5/2-j 

P  V 

y 

-(i 

f.) 

(11a) 


h"(T*) 


KU  (i.) 


(lib) 


gives 


J(y'.5,x) 


cm<t> 

U'aC4 


exp  h(T*)  ] 


00  1 

h  •*»  <- 


where 


_ 4 _ f/2 

7  £(T*)2J  u'aC4 

1/2 


exp  [-r2  h(T*) ] 


.  f*(x*) 


7^(T*)2 


y'a(T*)£(T*r 


r  5/ 2il 

exp 

2  .  *.  2 
-  5^(T  > 

I1  -y  j 

(12b) 


(12c) 


r'(T*) 


My 


1/2  1/7 

5<T)  ■  (A)  - 


(12d) 
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and 

f *(T*)  <■  n(T*)  (m/2’:kT*)"3/2  expK(T*)2]  .  (12e) 

In  going  from  (12b)  to  (12c)  I  have  used  the  fact  that  CM(T)ctT-3^  2 ,  and  that 
u'ci(T*)i(T*)6  -  1,  which  follows  from  (10).  Equation  (12d)  is  valid  provided 
T  <  T*  <  T^  or,  equivalently,  provided  1  <  ct£  ^  <  (T^/T)^^  and  u'  >  (aC^)-^. 

if 

It  is  easily  shown  that  when  T  <  the  integral  term  calculated  above 
dominates  the  upper  boundary  term  in  equation  (6b).  Thus, 


f(-',v,T) 


r 

1/2  -1 

r  5_f 

4- 

(-'^4)t,t*  f*(T*)  exp 

-  |  '(T*)2 

[‘  -  <i*)  ] 

_7£(T*)2_ 

(13a) 


which  is  valid  provided 


1  < 


and  (a£^)  *  <  u'  <  1 


(13b) 


where  T  (u'v)  is  defined  in  equation  (10). 

Equation  (I3a)  is  readily  interpreted:  the  distribution  function  at  u',v 
and  T  is  equal  to  a  local  Maxwellian  distribution  evaluated  at  v  and  at  a 
higher  temperature  T*(y'v),  times  the  probability  exp[-  yc(T*)2  (1  -  (ij*)^^}] 
that  a  (P^v)  electron  "emitted"  at  T*  will  survive  a  flight  to  T,  times  the 
interval  |  7^(T*)2  in  a^out  InT  ,  over  which  ,  v)  electrons  are 

effectively  emitted.  Note  that  since  T*  is  independent  of  T,  f  acquires  a 
dependence  on  T  only  through  the  1  -  (j*)^2  factor  appearing  in  the  argument 
of  the  exponential.  For  large  C  this  dependence  is  quite  weak,  as  may  be  seen 
from  Table  I,  where  values  of  and  £(T*)  are  listed  for  several  values  of 
their  arguments. 
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Since  T*  is  the  temperature  at  which  a  suprathermal  (u  ,v)  electron  (u<0) 

suffered  its  last  thermalizing  collision,  it  is  seen  from  Table  I  that 

f(u,v,T)  becomes  increasingly  nonlocal  in  character  as  z,  increases  and  as 

Note  also  that  the  quantity  t(T*)  has  several  meanings:  i)  it  is  the 

value  of  ^  with  which  C-.v)  electrons  leave  T*;  from  Table  I  is  it  seen  that 

2.25  <£(T*)  <  3.25.  Thus  it  is  the  near-thermal  part  of  the  distribution 

which  contributes  electrons  to  lower  temperature  regions.  ii)  The  quantity 

f  4  tt  ~|^  7  2  -  1.34/r(T*)  is  the  range  of  InT,  centered  on  InT*,  which 

L7  r(T*)2  J 

contributes  (_,v)  electrons  to  lower  temperatures.  iii)  The  quantity 

exp  [-^(T*)2  {1  -  (t£*)^2}]  is  the  probability  that  a  (p  ,v)  electron 

"emitted"  near  T*  will  survive  a  flight  to  T.  Surprisingly  perhaps,  this 

probability  decreases  as  increases  (because  T*  increases). 

An  approximate  expression  for  the  distribution  function  which  is  valid 
*  /TjA7/2 

for  T  >  Tfa  (i.e.,  for  a £”  >  \x  /  and  u<  0)  is  easily  obtained  from  (9a)  by 
expanding  the  quantity  h(T')  in  a  Taylor  series  about  T'  m  T^.  Adding  the 
upper  boundary  term  to  the  result  I  find  that 

t 

/  T  \7 / 2  (r  \772 

f(u,5,T)  -  Cm(T)  (exp[-C2  h(Th)]  -  ^ j  exp(-?2  h(T)])[l  +  o(^-|  ] 


IT  \7/2 

aC6  >  (~)  ;  -1  <  u  <-  (aS6)-1 


(14a) 


where,  from  (9b),  h(T)  *  1  and 


T  /t  \  f?T  \5/2  /  5/ 2l" 

-fe)  b 


(14b) 


The  inequality  in  (14b)  follows  from  the  conditions  T  <  T^  and  T  >  T^.  Note 
that  as  £  increases,  and  Z*  decrease  and,  in  consequence,  f  tends  to 


it  it 

f  (T^).  Note  also  that  f(Tjj)  *  f  (Tjj),  as  required  by  the  boundary  condition 
(5b). 

An  approximate  expression  for  the  distribution  of  upward  (u  >  0) 
propagating  suprathermal  electrons  can  be  obtained  by  expanding  the  argument 
of  Che  exponential  in  Che  integrand  (in  6a)  about  T'  m  T.  The  result  may  be 
written 


f( 


f*(T) 


U)^T)  «  -  +  Cm(Tc)  expt-C"  q(Tc 

1  +  (£  ) 


)] 


1  - 


(T/Tc) 


5/2 


,  A  ,T*v7/2 
1  +  K-f  ) 


(15n) 


0  <y  <  1 


where 


q(T') 


[( 


1  )5/2 
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and  now 


<JV/2  • 


(15b) 


(15c) 


The  salient  features  of  (15a)  are  the  following.  First,  if  eithery-^  0 
or  E,  — >  0,  then  f  — >  f*(T).  Second,  f(Tc)  ■  f*(Tc)  in  accord  with  (5a). 
Third,  for  T  >  Tc  and  T*  >  T,  the  second  term  in  (15a)  is  small  compared  to 
the  first.  It  therefore  follows  that 


1 

J  f(u,C,T)  du 

o 


lnCl+q^6) 

a?6 


f*(T) 


(16) 


so  that  there  are  significantly  fewer  suprathermal  electrons  propagating  up 
the  temperature  gradient  than  in  LTE. 


The  contribution  of  downward  propagating  electrons  to  the  angle-averaged 

6  7/2 

distribution  in  the  velocity  range  1  <  etc-  <  (T^/T)  can  be  found  by 
integrating  (13)  over  u  between  the  ranges  -1  <.  u  <  Using  (10)  and 

(13)  and  defining  U  (£,T)  to  be  (a£^)“*,  I  find  that 


1 

(?,T)  -  \  J  f(-u,C,T)  du 


(17a) 


~f —  cm(t)  j  exP  ^2y2(1  -  7  y5)l 


(17b) 


Expanding  the  argument  of  the  exponential  about  y  ■  u0^7  yields 

5/2) 


i  +  o(^6)-5/7 


J  1  <  a£6  <  ~ 


7/2 


(17c) 


where  f*  is  a  local  Maxwellian,?  (T*)  »  (7*)*^  and  f*  *  T*(v,y  *  -1);  that 
is. 


?  -  «.«6)2/7 


(17d) 


Note  that  T*  ■  f  (v)  only.  By  using  the  relation  ai^lx-T*  *  1,  equation 
(17d)  can  be  interpreted  in  a  way  similar  to  that  in  which  equation  (13)  was 
interpreted. 


'It* 
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On  comparing  (17c)  with  (16)  it  is  seen  that  upward  propagating  electrons 
make  a  negligible  contribution  to  the  angle-averaged  distribution  for  ot£6  » 
1.  Thus, 


Two  aspects  of  these  formulas  are  worth  mentioning.  F;-st,  since  both  a  and  T 
are  proportional  to  T-^2,  while  CM  is  proportional  to  T-^2,  it  is  seen  that 
f°  depends  on  T  only  through  the  (ct£6)  term  in  the  exponential.  The 

spatial  dependence  of  the  tail  is  therefore  quite  weak,  as  seen  earlier  in 
Figure  2.  Second,  at  low  temperatures  (large  a),  the  primary  velocity 
dependence  in  (18b)  is  through  the  £  '  term;  the  exponential  dependence  is 

relatively  weak.  For  example,  in  going  from  jr  *  3  to  C  ■  20  at  T  •  5  x  10^  K 
(a  -  9.23  x  10“3),  the  exponential  term  decreases  by  a  factor  of  1.9  x  102, 
while  the  £~31^7  term  decreases  by  a  factor  of  4.5  x  103.  In  low  temperature 
regions  the  angle-averaged  distribution  is  therefore  nearly  a  power  law,  as 
was  also  seen  earlier  in  Figures  1  and  2. 

Multiplying  equation  (13)  by  M and  integrating  from  p  ■  -1  to  y  ■  -(o£^)-* 


yields 


^U.T)  = 


-  !  /«- 


y>£>T)  dy  . 


(19a) 


-  -  3  f°  (£,T)  ;  1  <  o£ 


rf> 


(19b) 


where  f°(  ,T)  is  given  by  (18a)  or  (18b).  Thus  is  also  nearly  a  spatially 
independent,  v-^*^  power  law  in  the  low  temperature  regions  of  the 
atmosphere.  This  fact  accounts  for  the  slow  falloff  of  the  heat  flux 
integrand  shown  in  Figure  3. 

3 .  Generalizations 

The  above  results  are  easily  generalized  to  the  case  in  which  the 
background  temperature  profile  is  of  the  form 

Tr  -  Thr  +  (Thr  -  Tcr)  £  ;  r  >  1  (20a) 


n(z)T(z)  *  constant 


The  results  are: 


ar(T) 


T 

T 


/  \U2 

’(C.T)  -  (r-l)^j 


f*(T*) 

^(f*) 


” 

,  r-r 

exp 

”  r-T  (I*>2 

■4) 

_ 

1  +  o  (ar£6) 


(¥) 


(20b) 


(21a) 


(21b) 

(21c) 
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where  the  last  equation  is  valid  for 

i  <  cxrc6  <(r)  (21d> 

The  case  r  -  3  is  of  special  interest,  for  here  ctf  -  constant,  and  both  f°  and 
f ^  reduce  to  pure  v~^  power  laws. 

IV.  Summary 

In  the  preceding  I  have  shown  that  under  inhomogeneous  conditions  in 
which  only  the  tail  of  the  electron  velocity  distribution  departs 
significantly  from  a  local  Maxwellian  distribution,  as  is  the  case  in  the 
solar  transition  region,  a  linearized  BGK  collision  term  with  collision 
frequency  proportional  to  v  gives  results  for  the  tail  of  the  distribution 
which  are  in  adequate  agreement  with  results  obtained  from  the  Landau 

i 

equation.  Approximate  analytical  expressions  were  obtained  which  show  that 
the  zeroth  and  first  angular  moments  of  the  distribution  (and,  undoubtedly, 
higher  order  moments)  are  nearly  power  laws  over  a  wide  velocity  range  in  the 
lower  temperature  parts  of  the  atmosphere. 

The  linear  BGK  model  discussed  here  does  not  preserve  the  collisional 
invariants  and  hence  is  not  useful  for  determining  the  kinetic  temperature 
profile  corresponding  to  energy  transfer  solely  by  conduction.  Nevertheless, 
the  model  does  unambiguously  show  that  under  appropriate  conditions  (say,  ot(T) 

_  *5  C 

>  10  over  a  substantial  temperature  range)  the  energy  flux  into  the  lower 
temperature  regions  is  carried  by  suprathermal  electrons  with  velocities  much 
larger  than  local  thermal  velocities.  This  is  in  contrast  to  the  classical 
prediction  that  the  energy  flux  is  carried  by  electrons  with  velocities  two  to 
three  times  local  thermal  velocities. 

This  work  was  supported  in  part  by  NASA  grants  NGL  05-020-272  and  NAGW-92 
and  Office  of  Naval  Research  Contract  N00014-75-C-0673 . 
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Table  I 

Parameters  in  BGK  Solutions 


T  -  5  x 

104  K 

T  -  3  x 

105  K 

a  -  9.26 

x  10~3 

a  -  3.77 

x  10_J 

€ 

T*(u  ,£) 

T 

£(T*) 

2l(u,0 

T 

'  (T*) 

3 

-1.0 

1.73 

2.28 

1.33 

2.61 

3 

-0.5 

1.42 

2.51 

1.09 

2.87 

4 

-1.0 

2.82 

2.38 

2.19 

2.72 

4 

-0.5 

2.31 

2.63 

1.80 

2.98 

5 

-1.0 

4.14 

2.46 

3.20 

2.81 

5 

-0.5 

3.40 

2.71 

2.63 

3.08 

7 

-1.0 

7.37 

2.58 

5.71 

2.93 

7 

-0.5 

6.04 

2.85 

4.68 

3.24 

10 

-1.0 

13.6 

2.71 

— 

— 

10 

-0.5 

11.2 

2.99 

. 

— 

15 

15 


-1.0 

-0.5 


27.2 

22.3 


2.87 

3.17 


Figure 


Figure  2 


Figure  3 


Figure  4 


“1 


Figure  Captions 

1  Dimensionless,  angle-averaged  distibution  function  versus  £  *  v^vth" 
Solid  lines,  Landau  equation  results.  Broken  lines,  BGK  results. 
Dashed  line,  a  local  Maxwellian  distribution. 

Dimensional,  angle-averaged  distribution  function  f°  versus  velocity 
measured  in  units  of  the  thermal  velocity  at  the  upper  boundary. 
Solid  lines  are  BGK  results;  dashed  lines  are  local  Maxvellians .  . 

Absolute  value  of  the  heat  flux  integrand  versus  velocity  mesured  in 
units  of  thermal  velocity  at  the  upper  boundary.  Solid  lines,  full 
BGK  result.  Dashed  line,  classical  transport  theory  prediction  for 
BGK  equation. 

Heat  flux  calculated  from  linear  BGK  equation  (normalized  to  its 

l 

classical  value)  versus  temperature.  Results  for  several  pressures 
are  shown. 
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